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Abstract 

We present comprehensive molecular dynamics (MD) results for 
the kinetics of surface-directed spinodal decomposition (SDSD) and 
surface enrichment (SE) in binary mixtures at wetting surfaces. We 
study the surface morphology and the growth dynamics of the wetting 
and enrichment layers. The growth law for the thickness of these layers 
shows a crossover from a diffusive regime to a hydrodynamic regime. 
We provide phenomenological arguments to understand this crossover. 



A rich class of physical phenomena occurs when a homogeneous binary 
(AB) mixture is placed in contact with a surface (S). In many cases, the 
surface has a preferential attraction for one of the components of the mix- 
ture (say A). If the mixture is immiscible, the homogeneous bulk is unstable 
to phase separation and segregates into growing A-rich and .B-rich domains. 
The surface becomes the origin of surface- directed spinodal decomposition 
(SDSD) waves which propagate into the bulk HJElElIllISlElIZlISlini HOI HI] - 
The system evolves into a partially wet (PW) or completely wet (CW) equilib- 
rium morphology, depending upon the relative interaction strengths among 
A, B and S [121 D21 HH 03] • I n the PW morphology, the interface between 
the A-rich and B-rich domains meets the surface at a contact angle deter- 
mined by Young's condition [IB]. In the CW morphology, the A-rich phase 
expels the B-rich phase from the surface, and the AB interface is parallel 
to the surface. On the other hand, if the mixture is miscible, the emergent 
morphology consists of a thin surface enrichment (SE) layer followed by an 
extended depletion layer [171 HH [191 120]. The kinetics of SDSD and SE are 
of great technological and scientific importance, and find application in the 
fabrication of nanostructures, layered materials, composites, etc. 

Many experiments on SDSD and SE involve polymer or fluid mixtures, 
where hydrodynamic effects play an important role. For bulk phase-separation 
kinetics, it is well known that hydrodynamics has a drastic effect on the in- 
termediate and late stages [22] • The coarsening domains with size L(t) 
show a power-law growth, L(t) ~ t^, with the exponent changing from 
1/3 (diffusive regime) to 1 (viscous hydrodynamic regime) to 2/3 (inertial 
hydrodynamic regime) [231 [2U [22]. However, there is no analogous under- 
standing of the role of fluid velocity fields in SDSD and SE [26, 27]. In this 
paper, we present comprehensive molecular dynamics (MD) results for the 
kinetics of SDSD and SE. Our MD simulations clearly demonstrate a sharp 
crossover from a diffusive regime to a viscous hydrodynamic regime in both 
SDSD and SE. These results are of great relevance for experimentalists in 
this area, and provide a framework to understand the observation of diverse 
growth exponents in experiments. 

We consider a binary (AB) fluid mixture of point particles confined in a 
box of volume V = L x L x D. Periodic boundary conditions are applied in 
the x and y directions. An impenetrable surface is present at z — 0, which 
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gives rise to an integrated Lennard- Jones (LJ) potential (a = A,B): 
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where n is the fluid density, and a is the LJ diameter. In Eq. ([T]), e r and e a are 
the strengths of the repulsive and attractive parts of the surface potential. 
We set 5a = 1 and 5b = 0, i.e., A particles are attracted at large distances, 
whereas B particles are only repelled. Further, z' = z + a/2 so that the 
singularity of u w (z) occurs at z = —a/2 (outside the box). A similar surface 
is present at z = D with z' = D + a/2 — z, and 5a = 5b = 0, i.e., both A and 
B particles are repelled. The simulation box corresponds to a semi-infinite 
geometry [TJ: the generalization to any other geometry is straightforward. 
The particles in the system interact with LJ potentials (a = A,B): 



u(rij) = 4e Q/3 
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where = |r*j — fj\. We set the interaction parameters as €aa = £bb = 
2eab = £■ The bulk phase diagram for this potential is well known [28 | [29 | 130] . 
We use the truncated LJ potential with r c = 2.5a - this potential is shifted 
and force-corrected [31J. We consider a critical composition with Na = Nb = 
N/2. The system is a high-density incompressible fluid with n = N/V = 1. 
The particles have equal masses (rriA = tub = 1); and we set a = l,e = 
1, k B = 1 so that the MD time unit is t = ^Jma 2 / (48e) = l/v 7 ^- 

The MD runs were performed using the Verlet velocity algorithm [32] with 
a time-step At = 0.07 in MD units. We maintain the temperature (T) via the 
Nose- Hoover thermostat which is known to preserve hydrodynamics [321 123] . 
The homogeneous initial state for a run was prepared by equilibrating the 
system at high T with periodic boundary conditions in all directions. At 
time t — 0, the system is quenched to T < T c for SDSD (T c ~ 1.423) or 
T > T c for SE, and the surfaces are introduced at z = 0, D. We can also 
consider initial conditions where the mixture at high T has been equilibrated 
with the surfaces at t < 0. In that thin enrichment layer arises at 

z — - this layer vanishes at very high values of T. When the system is 
quenched, this very thin initial layer rapidly becomes part of the wetting or 
enrichment layers which grow from the surface. 

First, we present results for the kinetics of SDSD. The system size was 
L 2 xD with L = 48, D = 48 (N = 110592). The surface potential parameters 
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were e r = 0.5, e a = 0.6, and the quench temperature was T = 1.0 ~ 0.7T C . 
Equilibrium wetting phenomena for this model have been studied by Das 
and Binder [33J, who observed a first-order wetting transition at the temper- 
atures of interest. The above parameters correspond to a CW morphology in 
equilibrium. In Fig.[fl we show the laterally-averaged depth profiles [ip aY (z, t) 
vs. z] and the evolution snapshots (inset). The snapshot at t — 2800 shows 
the formation of an A-rich wetting layer at the surface (z = 0) in conjunction 
with phase separation in the bulk. This wetting layer propagates into the 
bulk, as seen from the depth profiles. The order parameter ip(r,t) is defined 
from the local densities as if)(r,t) = {tia — ?1b)/(^a + ns)- The quantity 
V'av(^j^) is obtained by averaging ip(r,t) in the x,y directions (parallel to 
the surface), and further averaging over 50 independent runs. In the bulk, 
the spinodal decomposition wave-vectors are randomly oriented - the above 
procedure yields ?/> av ~ 0. Near the surface, we see a structured morphology 
consisting of a wetting layer at the surface, depletion layer adjacent to it, 
etc. This layered structure propagates into the bulk. These SDSD profiles 
have been observed in many experiments on this problem [TJ [21 HU El [3] . 

An important characteristic of the SDSD profiles is the time-dependence 
of the wetting-layer thickness R\(t). This is defined as the first zero-crossing 
of the laterally-averaged depth profiles in Fig. [TJ In Fig. [21 we plot Ri(t) 
vs. t for the evolution shown in Fig. [TJ The growth dynamics is power-law, 
R^t) ~ t e , with 9 ~ 1/3 for t < 2000 and 9 ~ 1 for t > 2000. For these 
system sizes, we can go up to t ~ 3000 before encountering finite-size effects 
due to the lateral domain size becoming an appreciable fraction of the system 
size L. 

How can one understand this crossover? At early times, the wetting layer 
grows by the diffusive transport of A from bulk domains of size R ~ (at) 1 ^ 3 
(with chemical potential /i ~ cr/R, u being the surface tension) to the flat 
surface layer of size ~ oo (with \i ~ 0). If we neglect the very early potential- 
dependent growth regime [3l], we have 

dt Rh RRi 

where h ~ R± is the thickness of the depletion layer (see Fig. [TJ). From 
Eq. 0, we readily obtain the Lifshitz-Slyozov (LS) growth law: Ri ~ (crt) 1 ^ 3 . 
At later times, bulk tubes establish contact with the wetting layer and ma- 
terial is pumped hydrodynamically to the surface. The subsequent growth 
is analogous to that in phase separation of fluids - we expect Ri(t) ~ t 
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(viscous hydrodynamic regime) which crosses over to Ri(t) ~ t 2/3 (inertial 
hydrodynamic regime). The latter stage is presently not accessible via MD 
simulations [25j, due to computational limitations. However, our results for 
wetting-layer dynamics in Fig. [2] appear to access the viscous regime, albeit 
in a limited time-window. We remark that the crossover time (t c ~ 2000) is 
consistent with that reported by Ahmad et al. [25] in an MD simulation of 
bulk phase separation with similar parameter values. Clearly, we need sub- 
stantially larger system sizes to obtain an extended regime of linear growth 
in the viscous regime. Notice that the crossover in Fig. [2] is quite sharp, sug- 
gesting that there is a rapid pumping of material to the wetting layer when 
the bulk tubes first make contact. 

Second, we present results for the kinetics of SE. In this case, the system 
size was L 2 x D with L = 32, D = 64 (N = 66536). As the bulk remains 
homogeneous, the lateral size L (in the x, y directions) is not severely con- 
strained. However, in the direction perpendicular to the surface at z — 0, we 
need sufficiently large D to ensure decay of the enrichment profiles as z — > D. 
For the range of times studied here (t < 7000), test runs with other linear di- 
mensions showed that D = 64 is large enough to eliminate finite-size effects. 
In Fig. [3J we show the laterally-averaged profiles and an evolution snapshot 
(inset) for the kinetics of SE. The quench temperature was T = 2.0 ~ 1.41T C . 
The depth profiles were obtained using the same procedure as for SDSD, as 
an average over 50 independent runs. The surface potential parameters were 
e r = 0.5 and e a = 3.0. 

As expected, the morphology for SE is quite different from that in the 
SDSD case (cf. Fig. [1]). There is a thin enrichment layer of A at the surface. 
Due to the conservation of the order parameter, there must be a correspond- 
ing depletion layer which decays to ^ av ~ in the bulk. These profiles are in 
agreement with the experimental observations of Jones et al. [IB] on blends 
of deuterated and protonated polystyrene, and the experimental results of 
Mouritsen [35] on biopolymer mixtures. Notice that similar profiles are seen 
for SDSD if the system is quenched to the metastable region of the phase 
diagram [31]. The evolution dynamics in that case is analogous to the SE 
problem, as long as droplets are not nucleated in the system. 

For the case with diffusive dynamics, Binder and Frisch [36] and Frisch 
et al. [37] 138] have studied the morphology of SE profiles in the framework 
of a linear theory. They find that the SE profiles have a double-exponential 
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form: 



^(z,t)~B_(t) e~ z ^- {t) -B + (t) e" 2 ^ W , (4) 

with amplitudes B_(t), B + (t) > 0. The quantities B_(t) and £_(£) rapidly 
saturate to their equilibrium values a\ and &i, which depend on the surface 
potential. The other length scale grows diffusively with time, and B + (t) 
shows a corresponding decay: 

B + {t)^a 2 r^\ £ + (t)~M 1/2 . (5) 

The conservation constraint requires that = B + C, + . In Fig. [31 we see 

that the double-exponential function in Eq. (j3J) describes the SE profiles very 
well. 

What about the parameters of the double-exponential profile? As for the 
diffusive case, we find that -B_ and £_ rapidly saturate to their equilibrium 
values and may be treated as static quantities. The growing length scale 
is - we plot £ + vs. t for different surface field strengths in Fig. H](a). 

The early-time dynamics is consistent with diffusive growth (£ + ~ t 1 / 2 ), but 
again there is a crossover at t ~ t c to a hydrodynamic regime (£ + ~ t). 
The crossover is also seen in the thickness of the SE layer. For the profile 
in Eq. (jlj), the zero is located at Z (t) ~ ln(£>_/£> + ). Thus, we expect 
Z (t) ~ hit in both the time-regimes of Fig. Hta), but the slope should be 
steeper for t > t c . This is precisely the behavior seen in Fig. S](b), where 
we plot our MD results for Z (t) vs. t on a log-linear scale. This confirms 
that there is a crossover in the growth exponent from the diffusive regime 
(6 ~ 0.5) to the hydrodynamic regime (8 ~ 1.0). 

To understand the crossover in Fig. HI consider the dimensionless evolu- 
tion equation for the order parameter in the presence of a fluid velocity field 

^(f,t) = V 2 /^-v- VV, (6) 

where the chemical potential (for T > T c ) is /i = ip + ip 3 — (l/2)V 2 ip. For 
the SE problem, the system is homogeneous in the directions parallel to the 
surface. Thus, we set i/j(r,t) ~ ip(z,t) and v z (r, t) ~ v z (z,t). Then, Eq. flS]) 
becomes 
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We use the functional form of ip(z,t) in Eq. (J3J) to estimate the dominant 
contribution to the various terms in Eq. ([7]) at z ~ i.e., far from the 

surface [39J. We have 

dip 1 d£+ <9 2 /i 1 dip v z 

rsj rs_/ 7) ( N 1 

e+ dt ' ei' dz e+' 

For T > T c , the bulk is homogeneous and there is no large-scale structure 
formation in the composition or velocity fields. As the system is incom- 
pressible, the velocity field obeys V • v = or dv z /dz = in the laterally 
homogeneous case, so that v z ~ constant. Notice that v z < as there is 
a net current of the preferred component A towards the enrichment layer. 
This current is dissipated at the surface by the formation of locally inhomo- 
geneous structures. At early times, the diffusive term in Eq. (JTj) dominates, 
yielding £ + ~ i 1//2 . At late times, the convective term in Eq. ([7]) is dominant, 
giving a crossover to £ + ~ v z t. The precise dependence of the crossover time 
on various physical parameters can be estimated by considering the dimen- 
sional version of Eq. (JH]) in conjunction with the Navier-Stokes equation for 
the velocity field [211 [22]- 

In summary, we have presented comprehensive MD results for the kinetics 
of SDSD and SE in binary mixtures at wetting surfaces. Both the SDSD 
wetting layer and the SE layer show a crossover from a diffusive regime (with 
Ri ~ t 1 / 3 and £ + ~ t 1 ^ 2 ) to a hydrodynamic regime (with Ri ~ t and ~ t). 
These crossovers can be understood by simple phenomenological arguments. 
Our MD results are of great experimental relevance as most studies of these 
problems are done on polymer or fluid mixtures. We hope that our results in 
this paper will provoke fresh experimental interest in this problem, and our 
theoretical results will be subjected to experimental confirmation. 
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Figure 1: Kinetics of surface-directed spinodal decomposition (SDSD) in a 
binary (AB) Lennard- Jones (LJ) mixture. An impenetrable surface (located 
at z — 0) attracts the A-particles. The surface field strengths are e r = 
0.5, e a = 0.6, and the temperature is T = 1.0 ~ 0.7 T c . The other simulation 
details are given in the text. Main figure: Laterally- averaged order parameter 
profiles [ip av (z,t) vs. z\ at t = 350,700,2800 MD units. Inset: Evolution 
snapshots at t = 700, 2800. The A-particles are marked gray and the B- 
particles are marked black. 
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Figure 2: Time-dependence of the wetting-layer thickness [Ri(t) vs. t] of the 
SDSD profiles on a log- log scale. The straight lines have slope 1/3 (diffusive 
regime) and 1 (viscous hydrodynamic regime), respectively. 
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Figure 3: Kinetics of surface enrichment (SE) for an AB mixture in contact 
with a surface (located at z — 0) which attracts the A-particles. We set 
e r = 0.5, e a = 3.0, and T = 2.0 ~ 1.41 T c . The other simulation details are 
provided in the text. Main figure: Plot of ip w (z, t) vs. z at t = 280, 1400, 7000 
MD units. The double-exponential fits for the SE profiles are shown as solid 
lines. Inset: Evolution snapshot at t — 7000. The A-particles are marked 
gray and the 5-particles are marked black. 
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Figure 4: Evolution kinetics of SE profiles, (a) Time-dependence of £ + on 
a log-log scale for e r = 0.5 and e a = 1.2,2.0,3.0. The lines of slope 1/2 
and 1 correspond to diffusive and hydrodynamic growth, respectively, (b) 
Time-dependence of the zero-crossing Z (t) on a log-linear scale. 
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